Increased flexibility of the SARS-CoV-2 RNA-binding site causes resistance to remdesivir.

Mutations continue to accumulate within the SARS-CoV-2 genome, and the ongoing epidemic has shown no signs of ending. It is critical to predict problematic mutations that may arise in clinical environments and assess their properties in advance to quickly implement countermeasures against future variant infections. In this study, we identified mutations resistant to remdesivir, which is widely administered to SARS-CoV-2-infected patients, and discuss the cause of resistance. First, we simultaneously constructed eight recombinant viruses carrying the mutations detected in in vitro serial passages of SARS-CoV-2 in the presence of remdesivir. We confirmed that all the mutant viruses didn't gain the virus production efficiency without remdesivir treatment. Time course analyses of cellular virus infections showed significantly higher infectious titers and infection rates in mutant viruses than wild type virus under treatment with remdesivir. Next, we developed a mathematical model in consideration of the changing dynamic of cells infected with mutant viruses with distinct propagation properties and defined that mutations detected in in vitro passages canceled the antiviral activities of remdesivir without raising virus production capacity. Finally, molecular dynamics simulations of the NSP12 protein of SARS-CoV-2 revealed that the molecular vibration around the RNA-binding site was increased by the introduction of mutations on NSP12. Taken together, we identified multiple mutations that affected the flexibility of the RNA binding site and decreased the antiviral activity of remdesivir. Our new insights will contribute to developing further antiviral measures against SARS-CoV-2 infection.

Smoking enhances the expression of angiotensin-converting enzyme 2 involved in the efficiency of severe acute respiratory syndrome coronavirus 2 infection.

Smoking is one of the risk factors most closely related to the severity of coronavirus disease 2019 (COVID-19). However, the relationship between smoking history and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infectivity is unknown. In this study, we evaluated the ACE2 expression level in the lungs of current smokers, ex-smokers, and nonsmokers. The ACE2 expression level of ex-smokers who smoked cigarettes until recently (cessation period shorter than 6 months) was higher than that of nonsmokers and ex-smokers with a long history of nonsmoking (cessation period longer than 6 months). We also showed that the efficiency of SARS-CoV-2 infection was enhanced in a manner dependent on the angiotensin-converting enzyme 2 (ACE2) expression level. Using RNA-seq analysis on the lungs of smokers, we identified that the expression of inflammatory signaling genes was correlated with ACE2 expression. Notably, with increasing duration of smoking cessation among ex-smokers, not only ACE2 expression level but also the expression levels of inflammatory signaling genes decreased. These results indicated that smoking enhances the expression levels of ACE2 and inflammatory signaling genes. Our data suggest that the efficiency of SARS-CoV-2 infection is enhanced by smoking-mediated upregulation of ACE2 expression level.

Designing isolation guidelines for COVID-19 patients with rapid antigen tests.

Appropriate isolation guidelines for COVID-19 patients are warranted. Currently, isolating for fixed time is adopted in most countries. However, given the variability in viral dynamics between patients, some patients may no longer be infectious by the end of isolation, whereas others may still be infectious. Utilizing viral test results to determine isolation length would minimize both the risk of prematurely ending isolation of infectious patients and the unnecessary individual burden of redundant isolation of noninfectious patients. In this study, we develop a data-driven computational framework to compute the population-level risk and the burden of different isolation guidelines with rapid antigen tests (i.e., lateral flow tests). Here, we show that when the detection limit is higher than the infectiousness threshold values, additional consecutive negative results are needed to ascertain infectiousness status. Further, rapid antigen tests should be designed to have lower detection limits than infectiousness threshold values to minimize the length of prolonged isolation.

Estimation of timing of infection from longitudinal SARS-CoV-2 viral load data: mathematical modelling study

Multiple waves of the COVID-19 epidemic have hit most countries by the end of 2021. Most of those waves are caused by emergence and importation of new variants. To prevent importation of new variants, combination of border control and contact tracing is essential. However, the timing of infection inferred by interview is influenced by recall bias and hinders the contact tracing process. We propose a novel approach to infer the timing of infection, by employing a within-host model to capture viral load dynamics after the onset of symptoms. We applied this approach to ascertain secondary transmission which can trigger outbreaks. As a demonstration, the 12 initial reported cases in Singapore, which were considered as imported because of their recent travel history to Wuhan, were analyzed to assess whether they are truly imported. Our approach suggested that 6 cases were infected prior to the arrival in Singapore, whereas other 6 cases might have been secondary local infection. Three among the 6 potential secondary transmission cases revealed that they had contact history to previously confirmed cases. Contact trace combined with our approach using viral load data could be the key to mitigate the risk of importation of new variants by identifying cases as early as possible and inferring the timing of infection with high accuracy.

Different efficacies of neutralizing antibodies and antiviral drugs on SARS-CoV-2 Omicron subvariants, BA.1 and BA.2.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron subvariant BA.2 has spread in many countries, replacing the earlier Omicron subvariant BA.1 and other variants. Here, using a cell culture infection assay, we quantified the intrinsic sensitivity of BA.2 and BA.1 compared with other variants of concern, Alpha, Gamma, and Delta, to five approved-neutralizing antibodies and antiviral drugs. Our assay revealed the diverse sensitivities of these variants to antibodies, including the loss of response of both BA.1 and BA.2 to casirivimab and of BA.1 to imdevimab. In contrast, EIDD-1931 and nirmatrelvir showed a more conserved activities to these variants. The viral response profile combined with mathematical analysis estimated differences in antiviral effects among variants in the clinical concentrations. These analyses provide essential evidence that gives insight into variant emergence's impact on choosing optimal drug treatment.

Incomplete antiviral treatment may induce longer durations of viral shedding during SARS-CoV-2 infection.

The duration of viral shedding is determined by a balance between de novo infection and removal of infected cells. That is, if infection is completely blocked with antiviral drugs (100% inhibition), the duration of viral shedding is minimal and is determined by the length of virus production. However, some mathematical models predict that if infected individuals are treated with antiviral drugs with efficacy below 100%, viral shedding may last longer than without treatment because further de novo infections are driven by entry of the virus into partially protected, uninfected cells at a slower rate. Using a simple mathematical model, we quantified SARS-CoV-2 infection dynamics in non-human primates and characterized the kinetics of viral shedding. We counterintuitively found that treatments initiated early, such as 0.5 d after virus inoculation, with intermediate to relatively high efficacy (30-70% inhibition of virus replication) yield a prolonged duration of viral shedding (by about 6.0 d) compared with no treatment.

Revisiting the guidelines for ending isolation for COVID-19 patients.

Since the start of the COVID-19 pandemic, two mainstream guidelines for defining when to end the isolation of SARS-CoV-2-infected individuals have been in use: the one-size-fits-all approach (i.e. patients are isolated for a fixed number of days) and the personalized approach (i.e. based on repeated testing of isolated patients). We use a mathematical framework to model within-host viral dynamics and test different criteria for ending isolation. By considering a fixed time of 10 days since symptom onset as the criterion for ending isolation, we estimated that the risk of releasing an individual who is still infectious is low (0-6.6%). However, this policy entails lengthy unnecessary isolations (4.8-8.3 days). In contrast, by using a personalized strategy, similar low risks can be reached with shorter prolonged isolations. The obtained findings provide a scientific rationale for policies on ending the isolation of SARS-CoV-2-infected individuals.

Detection of significant antiviral drug effects on COVID-19 with reasonable sample sizes in randomized controlled trials: A modeling study.

BACKGROUND: Development of an effective antiviral drug for Coronavirus Disease 2019 (COVID-19) is a global health priority. Although several candidate drugs have been identified through in vitro and in vivo models, consistent and compelling evidence from clinical studies is limited. The lack of evidence from clinical trials may stem in part from the imperfect design of the trials. We investigated how clinical trials for antivirals need to be designed, especially focusing on the sample size in randomized controlled trials. METHODS AND FINDINGS: A modeling study was conducted to help understand the reasons behind inconsistent clinical trial findings and to design better clinical trials. We first analyzed longitudinal viral load data for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) without antiviral treatment by use of a within-host virus dynamics model. The fitted viral load was categorized into 3 different groups by a clustering approach. Comparison of the estimated parameters showed that the 3 distinct groups were characterized by different virus decay rates (p-value < 0.001). The mean decay rates were 1.17 d-1 (95% CI: 1.06 to 1.27 d-1), 0.777 d-1 (0.716 to 0.838 d-1), and 0.450 d-1 (0.378 to 0.522 d-1) for the 3 groups, respectively. Such heterogeneity in virus dynamics could be a confounding variable if it is associated with treatment allocation in compassionate use programs (i.e., observational studies). Subsequently, we mimicked randomized controlled trials of antivirals by simulation. An antiviral effect causing a 95% to 99% reduction in viral replication was added to the model. To be realistic, we assumed that randomization and treatment are initiated with some time lag after symptom onset. Using the duration of virus shedding as an outcome, the sample size to detect a statistically significant mean difference between the treatment and placebo groups (1:1 allocation) was 13,603 and 11,670 (when the antiviral effect was 95% and 99%, respectively) per group if all patients are enrolled regardless of timing of randomization. The sample size was reduced to 584 and 458 (when the antiviral effect was 95% and 99%, respectively) if only patients who are treated within 1 day of symptom onset are enrolled. We confirmed the sample size was similarly reduced when using cumulative viral load in log scale as an outcome. We used a conventional virus dynamics model, which may not fully reflect the detailed mechanisms of viral dynamics of SARS-CoV-2. The model needs to be calibrated in terms of both parameter settings and model structure, which would yield more reliable sample size calculation. CONCLUSIONS: In this study, we found that estimated association in observational studies can be biased due to large heterogeneity in viral dynamics among infected individuals, and statistically significant effect in randomized controlled trials may be difficult to be detected due to small sample size. The sample size can be dramatically reduced by recruiting patients immediately after developing symptoms. We believe this is the first study investigated the study design of clinical trials for antiviral treatment using the viral dynamics model.

Time variation in the probability of failing to detect a case of polymerase chain reaction testing for SARS-CoV-2 as estimated from a viral dynamics model.

Viral tests including polymerase chain reaction (PCR) tests are recommended to diagnose COVID-19 infection during the acute phase of infection. A test should have high sensitivity; however, the sensitivity of the PCR test is highly influenced by viral load, which changes over time. Because it is difficult to collect data before the onset of symptoms, the current literature on the sensitivity of the PCR test before symptom onset is limited. In this study, we used a viral dynamics model to track the probability of failing to detect a case of PCR testing over time, including the presymptomatic period. The model was parametrized by using longitudinal viral load data collected from 30 hospitalized patients. The probability of failing to detect a case decreased toward symptom onset, and the lowest probability was observed 2 days after symptom onset and increased afterwards. The probability on the day of symptom onset was 1.0% (95% CI: 0.5 to 1.9) and that 2 days before symptom onset was 60.2% (95% CI: 57.1 to 63.2). Our study suggests that the diagnosis of COVID-19 by PCR testing should be done carefully, especially when the test is performed before or way after symptom onset. Further study is needed of patient groups with potentially different viral dynamics, such as asymptomatic cases.

Potential anti-COVID-19 agents, cepharanthine and nelfinavir, and their usage for combination treatment.

Antiviral treatments targeting the coronavirus disease 2019 are urgently required. We screened a panel of already approved drugs in a cell culture model of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and identified two new agents having higher antiviral potentials than the drug candidates such as remdesivir and chroloquine in VeroE6/TMPRSS2 cells: the anti-inflammatory drug cepharanthine and human immunodeficiency virus protease inhibitor nelfinavir. Cepharanthine inhibited SARS-CoV-2 entry through the blocking of viral binding to target cells, while nelfinavir suppressed viral replication partly by protease inhibition. Consistent with their different modes of action, synergistic effect of this combined treatment to limit SARS-CoV-2 proliferation was highlighted. Mathematical modeling in vitro antiviral activity coupled with the calculated total drug concentrations in the lung predicts that nelfinavir will shorten the period until viral clearance by 4.9 days and the combining cepharanthine/nelfinavir enhanced their predicted efficacy. These results warrant further evaluation of the potential anti-SARS-CoV-2 activity of cepharanthine and nelfinavir.

A quantitative model used to compare within-host SARS-CoV-2, MERS-CoV, and SARS-CoV dynamics provides insights into the pathogenesis and treatment of SARS-CoV-2.

The scientific community is focused on developing antiviral therapies to mitigate the impacts of the ongoing novel coronavirus disease 2019 (COVID-19) outbreak. This will be facilitated by improved understanding of viral dynamics within infected hosts. Here, using a mathematical model in combination with published viral load data, we compare within-host viral dynamics of SARS-CoV-2 with analogous dynamics of MERS-CoV and SARS-CoV. Our quantitative analyses using a mathematical model revealed that the within-host reproduction number at symptom onset of SARS-CoV-2 was statistically significantly larger than that of MERS-CoV and similar to that of SARS-CoV. In addition, the time from symptom onset to the viral load peak for SARS-CoV-2 infection was shorter than those of MERS-CoV and SARS-CoV. These findings suggest the difficulty of controlling SARS-CoV-2 infection by antivirals. We further used the viral dynamics model to predict the efficacy of potential antiviral drugs that have different modes of action. The efficacy was measured by the reduction in the viral load area under the curve (AUC). Our results indicate that therapies that block de novo infection or virus production are likely to be effective if and only if initiated before the viral load peak (which appears 2-3 days after symptom onset), but therapies that promote cytotoxicity of infected cells are likely to have effects with less sensitivity to the timing of treatment initiation. Furthermore, combining a therapy that promotes cytotoxicity and one that blocks de novo infection or virus production synergistically reduces the AUC with early treatment. Our unique modeling approach provides insights into the pathogenesis of SARS-CoV-2 and may be useful for development of antiviral therapies.

Estimation of the incubation period of COVID-19 using viral load data.

The incubation period, or the time from infection to symptom onset, of COVID-19 has usually been estimated by using data collected through interviews with cases and their contacts. However, this estimation is influenced by uncertainty in the cases' recall of exposure time. We propose a novel method that uses viral load data collected over time since hospitalization, hindcasting the timing of infection with a mathematical model for viral dynamics. As an example, we used reported data on viral load for 30 hospitalized patients from multiple countries (Singapore, China, Germany, and Korea) and estimated the incubation period. The median, 2.5, and 97.5 percentiles of the incubation period were 5.85 days (95 % CI: 5.05, 6.77), 2.65 days (2.04, 3.41), and 12.99 days (9.98, 16.79), respectively, which are comparable to the values estimated in previous studies. Using viral load to estimate the incubation period might be a useful approach, especially when it is impractical to directly observe the infection event.